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This paper continues a numerical investigation of the statistical properties of orbits evolved in 
'frozen,' time- independent A^-body realisations of smooth time- independent density distributions p 
corresponding to both integrable and nonintegrable potentials, allowing for 10^'^ ^ N < 10*'^. The 
principal focus is on distinguishing between, and quantifying, the effects of graininess on initial 
conditions corresponding, in the continuum limit, to regular and chaotic orbits. Ordinary Lya- 
punov exponents x do not provide a useful diagnostic for distinguishing between regular and chaotic 
behaviour. Frozen- orbits corresponding in the continuum limit to both regular and chaotic charac- 
teristics have large positive x even though, for large A'^, the 'regular' frozen- A'^ orbits closely resemble 
regular characteristics in the smooth potential. Alternatively, viewed macroscopically both 'regular' 
and 'chaotic' frozen-A'^ orbits diverge eis a power law in time from smooth orbits with the same 
initial condition. There does, however, exist an important difference between 'regular' and 'chaotic' 
frozen-A'' orbits: For regular orbits, the time scale associated with this divergence 
with to a characteristic dynamical, or crossing time; for chaotic orbits to ~ {In N)tD- Convergence 
towards the continuum limit is much slower for chaotic orbits. At least for A'' > 10"^ or so, clear 
distinctions exist between phase mixing of initially localised orbit ensembles which, in the continuum 
limit, exhibit regular versus chaotic behaviour. Regular ensembles evolved in a frozen-A' density dis- 
tribution diverge as a power law in time, albeit more rapidly than ensembles evolved in the smooth 
distribution. Alternatively, chaotic ensembles diverge in a fashion that is roughly exponential, albeit 
at a larger rate than that associated with the exponential divergence of the corresponding ensemble 
evolved in the smooth p. For both regular and chaotic ensembles, finite- A' effects are well mimicked, 
both qualitatively and quantitatively, by energy-conserving white noise with amplitude rj cx 1/N. 
This suggests strongly that earlier investigations of the effects of low amplitude noise on phase space 
transport in smooth potentials are directly relevant to real physical systems. 



PACS number(s): 05.60.-hw, 51.10.+y, 05.40.+j 
I. MOTIVATIONS AND EXPECTATIONS 

This is the second in a series of papers, the aim of 
which is to understand the role of discreteness effects, 
i.e., graininess, in the gravitational A'^-body problem. 
Particular emphasis is given to the meaning of chaos and 
various manifestations of chaotic behaviour. As discussed 
in the first paper Q (hereafter Paper I) , this problem can 
be divided into two separate components, namely first 
understanding how graininess alters the motions of rep- 
resentative orbits in a fixed gravitational potential and 
only then considering how these changes are manifested 
in the context of a fully self-consistent A^-body evolution. 
As in Paper I, the focus here is on the former issue. 

One is thus led naturally to effect a statistical compari- 
son between (1) orbits evolved in a frozen iV-body density 
distribution generated by randomly sampling some spec- 
ified smooth density distribution p and (2) orbits evolved 
in the smooth potential $ related to p by Poisson's equa- 
tion, V^* = AnGp. 

In this setting, the present paper has two specific ob- 
jectives, namely (1) to implement precise, quantitative 



distinctions between 'regular' and 'chaotic' behaviour in 
such frozen-iV systems; and (2) to determine the extent 
to which discreteness effects can be mimicked successfully 
by a suitably defined 'noise' acting on orbits evolving in 
an otherwise smooth potential. 

Implementing useful distinctions between regular and 
chaotic behaviour is not completely trivial. For example, 
ordinary Lyapunov exponents computed for individual 
orbits do not provide a useful diagnostic. Even for den- 
sity distributions corresponding to integrable potentials, 
A'^-body orbits have large positive Lyapunov exponents. 
Moreover, even though there is a precise sense in which, 
as N increases, frozen- A^ orbits come to more closely re- 
semble orbits in the smooth potential, the values of these 
exponents do not decrease systematically with increasing 
N§§. 

Viewed macroscopically, a frozen- A^ orbit and a 
smooth orbit evolved from the same initial condition in 
density distributions corresponding to an integrable po- 
tential will typically diverge linearly in time on a time 
scale to fx A^^/^. As discussed in Paper I, this superfi- 
cially surprising result would appear to reflect the fact 
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that the chaos is associated with a large number of en- 
counters with neighbouring particles, each of very short 
duration, which tend to cancel systematically so as to 
have a comparatively minor macroscopic effect. 

For frozen- systems with very small N, one would an- 
ticipate that close encounters are relatively stronger and 
'more macroscopic' in nature than is the case for larger 
N, so that clean macroscopic distinctions between regu- 
lar and chaotic behaviour might be difficult to implement. 
However, for larger TV such close encounters become pro- 
gressively 'more microscopic', so that one might expect 
comparatively clear cut distinctions to exist. 

Following the pioneering work of Chandrasekhar in the 
1940's 1^ one might expect that discreteness effects 
will act in much the same way as friction and white noise, 
so that they can be modelled in the context of a Fokker- 
Planck description. To the extent that this intuition is 
correct, earlier work probing the effects of friction and 
noise on smooth potential orbits translate into precise 
predictions as to the expected effects of graininess. 

The amplitude rj associated with friction and noise de- 
fines a characteristic relaxation time tj^ ~ l/rj on which, 
e.g., the perturbation will induce significant changes in 
conserved quantities like energy. However, modelling dis- 
creteness effects as a sequence of close binary encounters 
which result in friction and noise yields the concrete pre- 
diction 1^ that til oc {N/ hiN)t]j, where is a charac- 
teristic dynamical, or crossing, time. One might, there- 
fore, anticipate that discreteness effects associated with a 
system comprised of N bodies can be reproduced by fric- 
tion and noise with characteristic amplitude 77 cx In N/N. 

It is well known that, when subjected to friction and 
white noise, orbits in an integrable potential and reg- 
ular orbits in a generic potential typically diverge as a 
power law from unperturbed orbits with the same initial 
condition. One might, therefore, expect that discrete- 
ness effects also induce a power law divergence between 
frozen- iV orbits and smooth characteristics with the same 
initial condition, and that the divergence time scale tc 
associated with a system of N bodies can be mimicked 
by noise with amplitude rj oc \nN/N. As noted in Paper 
I, frozen- iV orbits and smooth characteristics do indeed 
tend to diverge linearly on a time scale to oc N^/^to; 
and, as will be seen below, this linear divergence is well 
reproduced by an appropriately defined white noise with 
amplitude rj oc (Given the limited dynamical range 

in particle number for the simulations described in this 
paper - 10^'^ < N < 10^-^ - it would seem impossible 
to distinguish unambiguously between scalings 77 oc 1 /N 
and 770c In N/N. The simulations are consistent with 
both.) 

When subjected to friction and white noise, chaotic 
orbits behave very differently. Comparatively weak per- 
turbations typically induce an initial exponential diver- 
gence from the unperturbed orbit at a rate A that is com- 
parable to the largest (short time) Lyapunov exponent 
Xs for the unperturbed orbits [^. For stronger pertur- 
bations, the separation between perturbed and unper- 



turbed orbits quickly becomes 'macroscopic,' at which 
point the divergence become slower than exponential, al- 
beit still more rapid than what is observed for regular or- 
bits |6|. One might, therefore, expect that, when acting 
on chaotic initial conditions, discreteness effects associ- 
ated with a finite- iV system would induce (1) an initial 
exponential divergence at a rate A that is comparable to 
XSi a number typically much smaller than the 'true' Lya- 
punov exponent x associated with orbits in the frozen- TV 
system, followed by (2) a slower subexponential diver- 
gence which is still faster than the divergence associated 
with regular orbits. 

The simulations summarised in this paper provided 
unambiguous confirmation of the second of these ex- 
pectations. Viewed macroscopically, frozen-A^ orbits 
corresponding in the continuum limit to chaotic orbits 
typically diverge linearly from smooth orbits with the 
same initial condition; but for large N the time scale 
to oc {hiN)t]j associated with this divergence is much 
shorter than the time scale to oc N^^^tu associated with 
regular orbits. 

Unperturbed orbits evolved in a smooth integrable po- 
tential are multiply periodic and, as such, are charac- 
terised by Fourier spectra in which power is concen- 
trated at a countable set of discrete frequencies. Fric- 
tion and noise destroy this exact periodicity, resulting in 
a more complex Fourier spectrum. To the extent that 
the friction and noise are weak, the orbit should remain 
nearly regular and the spectrum should remain 'similar 
to' the spectrum associated with the unperturbed or- 
bit. However, when the friction and noise become larger 
in amplitude the orbit should become 'less nearly peri- 
odic,' and the spectrum should become 'more complex' 
than the spectrum associated with the unperturbed or- 
bit. In the same sense, one might expect that, as N 
decreases, frozen- orbits corresponding to regular or- 
bits in a smooth potential will become 'less regular' and 
be characterised by Fourier spectra that are 'more com- 
plex'. As described in Paper I, this intuition can be made 
precise by computing from an orbital time series such 
quantifiable measures of orbital complexity as (1) the 
number of frequencies in a discrete Fourier series that 
contain more than some fixed fraction j of the power in 
the peak frequency or (2) the minimum number of fre- 
quencies required to capture a fixed fraction k of the total 
power. 

By contrast, chaotic orbits in a smooth potential are 
in general aperiodic This implies that, even in a dis- 
crete time series, their power should be spread over a 
larger number of frequencies, so that the orbit will have 
substantially larger 'complexity' 0. For very small TV, 
where the qualitative character of the orbits is dominated 
by close encounters and distinctions between chaotic and 
regular motions are difficult to identify, one might ex- 
pect that frozen- orbits corresponding in the continuum 
limit to regular and chaotic orbits would have comparable 
complexities, much larger than the typical complexity as- 
sociated with a smooth regular orbit and larger even than 
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the complexity associated with a smooth chaotic orbit. 
As N increases, discreteness effects will presumably be- 
come less important and the complexity of both the regu- 
lar and chaotic frozen- orbits will decrease. For the case 
of regular frozen- iV orbits, the complexity should even- 
tually converge towards the comparatively small value 
associated with a smooth regular orbit. For the case of 
chaotic frozen-A orbits, the complexity should instead 
converge towards the substantially larger value associ- 
ated with a smooth chaotic orbit. 

To the extent that discreteness effects associated with 
a fixed number A can be successfully modelled in terms 
of a suitably defined noise with amplitude 77, one might 
also expect that the typical complexity n{N) associated 
with a frozen-A orbit with given A will be comparable 
to the complexity n{ri) associated with a smooth orbit 
with the same initial condition evolved in the presence 
of perturbations of amplitude 77 cx 1/A^. As described 
below, both these intuitions were in fact confirmed. 

Discreteness effects can also be explored in the context 
of the phase mixing of initially localised orbit ensembles, 
a subject which has received considerable attention in 
both galactic astronomy |^ and accelerator dy- 

namics [|l2| since Merritt and Valluri |l4| coined the 
term 'chaotic mixing' to characterise the much more effi- 
cient phase mixing associated with ensembles of chaotic 
orbits. 

Localised ensembles of regular initial conditions 
evolved in a smooth potential will initially diverge as a 
power law in time and, when viewed over much longer 
time scales, exhibit a coarse-grained evolution, again 
proceeding as a power law in time, towards a time- 
independent equilibrium state. The introduction of fric- 
tion and noise accelerates the original divergence, but 
that divergence still proceeds as a power law in time. 
Quantifying the later time evolution is more subtle be- 
cause the perturbations allow the orbits to access phase 
space regions which would otherwise be inaccessible. 
However, what is clear is that, as probed by various 
lower order moments, the ensemble evolves exponentially 
in time towards a 'well-mixed' state that manifests the 
symmetries of the unperturbed potential If, e.g., the 
potential admits a reflection symmetry, $(— x) = ^{x), 
the mean value (x) associated with the ensemble will con- 
verge exponentially towards zero. 

By contrast, localised ensembles of chaotic orbits 
evolved in a smooth potential initially diverge exponen- 
tially at a rate A that is comparable to the value of 
the largest Lyapunov exponent xs and, when viewed 
over somewhat longer time scales, exhibit a coarse- 
grained evolution, exponential in time, towards a time- 
independent, or nearly time-independent state. The rate 
A associated with this subsequent evolution is typically 
much smaller than A and is not directly related to xs, 
although loose correlations between A and xs often exist 
|6). Subjecting these same ensembles to friction and 
noise typically increases the rate A associated with the 
initial divergence, making the orbits behave in a fashion 



that is 'even more chaotic' Over sufficiently long time 
scales, these perturbations will eventually drive the en- 
semble towards a thermal state with a temperature Q set 
by the friction and noise. However, on time scales much 
shorter than the natural time scale tfi associated with the 
friction and noise, the ensemble will again evolve towards 
a nearly time-independent distribution; and, in many 
cases, the perturbations will increase the rate A asso- 
ciated with this convergence towards a near-equilibrium 

i- 

The simulations described in this paper demonstrate 
that, both qualitatively and quantitatively, discreteness 
effects impact phase mixing for both regular and chaotic 
orbits in much the same way as these perturbations. 

Perhaps the most important conclusion of the work de- 
scribed here is the apparent need to distinguish between 
two distinct 'types' of chaos. Short-range microscopic 
chaos associated with close encounters between individ- 
ual masses, is a ubiquitous phenomenon for the A^-body 
problem, which appears to be present irrespective of the 
bulk properties of the density distribution. In addition, 
however, there is the possibility of macroscopic chaos, 
easily identified in the continuum limit, which has spe- 
cific predictable implications for motions in the A^-body 
problem. 

Section II describes the potentials that were consid- 
ered and the numerical experiments that were performed. 
Section HI described the results derived for individual 
frozen- A?^ trajectories, determining how these results scale 
with N and demonstrating the extent to which they can 
be mimicked by a suitably defined white noise. Section 
IV then describes the results derived for the phase mix- 
ing of orbit ensembles, again considering both how things 
scale with A^ and the degree to which discreteness ef- 
fects can be mimicked by noise. Section V focuses on the 
possibility of transitions between 'regular' and 'chaotic' 
behaviour, a phenomenon which, especially for small A^, 
can be important for systems where the smooth poten- 
tial allows a coexistence of large measures of both regular 
and chaotic orbits. Section VI concludes by summaris- 
ing the principal conclusions and speculating on potential 
implications. 

II. DESCRIPTION OF THE NUMERICAL 
EXPERIMENTS 

The numerical experiments described in this paper fo- 
cused on the behaviour of orbits and orbit ensembles 
evolved in frozen-A^ realisations of four different time- 
independent density distributions. In the continuum 
limit, two of these correspond to integrable potentials; 
the other two correspond to potentials which admit large 
measures of chaos. These orbits and orbit ensembles were 
also compared with orbits with the same initial condi- 
tions evolved in the smooth potentials associated with 
the smooth density distributions, both with and without 
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noise. Each system was normalised to have mass M = 1, 
and units were so chosen that the gravitational constant 
G = 1. The four density distributions were the following: 
1. A spherically symmetric Plummer sphere, for which 
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This corresponds via Poisson's equation to a potential of 
the form 



$p(r) = - 



GM 
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This integrable system was considered extensively in Pa- 
per I. As in that paper, units were chosen such that b — I. 
2. A constant density triaxial ellipsoid, for which 
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For m < 1, this yields a potential of the form 
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where the frequencies Wa, t^b, are related to the axis 
values a, 6, c, in terms of incomplete elliptic integrals 
[ p5[ . Attention focused primarily on the specific param- 
eter values a — 1.95, 6 = 1.50, and c = 1.05, which 
imply $0 ~ - 2.91758, uja ~ 0.9759, « 0.8559, and 
ujc « 0.7161. 

3. A constant density ellipsoid perturbed by a supermas- 
sive black hole, this corresponding to a potential of the 
form 
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with e = 10 Attention here focused on a black hole 
mass Mbh = IQ-^ '^M « 0.0316228M, which yields ||l|] 
a potential for which, for orbits restricted energetically 
to m < 1, the phase space is largely chaotic. 
4. A triaxial generalisation of the Dehnen [|l^ potential, 
for which 



PD{r) 



M(3 - 7) 
47ra6c 



TO (1 + to) 



~(4-7) 



(2.7) 



Attention here focused on the parameter values j — 1 
and a ~ 1.0, 6 = -y/S. 0/8.0, and c = 0.5, values first 
considered by Merritt and Fridman ||l^ as a prototype for 
a cuspy triaxial galaxy. The phase space associated with 
this potential has been studied extensively |18| |J9| and 
it is known that, especially for low energies, there exist 
large measures of chaotic orbits. The smooth potential 



$£) associated with po cannot be written analytically, 
so that integrations in $£) were performed using a 32nd 
order Gauss- Legendre integration scheme written by C. 
Siopis as part of his Ph. D. dissertation l20[| . 

Independent of energy, the characteristical orbital time 
scale for motion in the constant density ellipsoid, with or 
without a black hole, corresponds to a dynamical time 
to ''^ 10. (The quantity 1/^/Gp ~ 3.6.) The experiments 
with the Plummer potential described in this paper were 
performed for intermediate energies for which, again, 
tj) ^ 10. The Dehnen potential exhibits a much larger 
degree of central concentration and, as such, to exhibits 
greater variability. 

The nonintegrable potential (2.6) is comparatively sim- 
ple in the sense that, for the energies and black hole 
masses considered here, almost all smooth orbits are 
chaotic. This implies that, when considering A^-body re- 
alisations, one need not be much concerned about the 
possibility of graininess converting an initially 'chaotic' 
orbit into a 'regular' orbit. By contrast, for the compar- 
atively low energies considered here the smooth Dehnen 
potential admits large measures of both regular and 
chaotic orbits, so that graininess could well induce nu- 
merous transitions between 'regular' and 'chaotic' be- 
haviour. 

Frozen- A^ density distributions of the form 
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were generated by randomly sampling the smooth density 
distributions p. These correspond to A'^-body potentials 
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which incorporate a tiny 'softening parameter' with value 
e < 10-3. 

Orbits were integrated in frozen-A^ realisations with 
102-5 < AT < 105-5. The integrations were performed with 
a variable time step integration scheme which was guar- 
anteed to conserve the energy of each particle to at least 
one part in 10"^. The energy of a typical orbit was 
conserved to within a few parts in 10^. Estimates of 
the largest (short time) Lyapunov exponents orbits were 
obtained in the usual way by simultaneously tracking 
the evolution of a small initial perturbation, periodically 
renormalised at fixed intervals At j2^. The short time 
Lyapunov exponents derived in this fashion typically ex- 
hibited rapid convergence towards near-constant values 
- much more rapid convergence than has been observed 
|]l6| ||l9[| for the smooth potentials - so that a time ~ 100 
dynamical, or crossing, times tu was sufficient to yield 
reasonable estimates. 

Some of the orbital data were Fourier analysed to de- 
termine their 'orbital complexity', defined as in |Q| as the 
number of frequencies required to capture a fixed frac- 
tion fc of the total power. More precisely, this entailed 
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determining for each orbit the quantities Ux, riy, and n^, 
defined, respectively, as the minimum number of frequen- 
cies required to capture a fixed fraction k of the power 
in each direction, and then assigning a total complexity 

n^nx + ny + Uz- (2-10) 

In order to obtain a reasonably sharp Fourier spec- 
trum, orbital data were typically recorded at intervals 
St = O.Oltu or less, and each orbit was represented by a 
time series containing at least 4096 points. 

Phase mixing in frozen-A'^ systems was explored by 
performing integrations of orbit ensembles comprised 
(mostly) of 1600 initial conditions localised within a 
phase space hypercube of size ~ lO^'^ the size of the ac- 
cessible phase space region. For each cell of initial con- 
ditions, experiments were typically repeated for several 
different frozen- iV density distributions, each sampling 
the same smooth p, so as to extract improved statistics. 
The experiments with N = 10^'^, 10^ and 10^'^ were each 
performed for six different frozen- iV density distributions, 
those with N = and 10*-^ for four distributions, and 
those for N = 10^ for two distributions. Because of com- 
putational constraints - each orbit with A'^ = 10^-^ aver- 
aged roughly three hours on a Pentium 200 workstation 
-, the long time ensemble integrations for N = 10^-^ 
were performed for only one ensemble, comprised of 800, 
rather than 1600, orbits. However, the first fifth of each 
integration, in many respects the most interesting, was 
repeated for a second frozen- iV distribution. 

To test the physical intuition that discreteness effects 
can be mimicked by friction and white noise in the con- 
text of a Fokker-Planck description, orbit ensembles were 
also evolved in the smooth potential in the presence of 
a suitably defined 'white noise.' Ordinarily, discreteness 
effects are modeled by considering a Langevin equation 
of the form || 

T dv 
-^^-\/,<i>-rj^+Fa, a = x,y,z. (2.11) 

Here rjdra/dt represents a dynamical friction and Fa rep- 
resents Gaussian white noise, which is characterised com- 
pletely by its first two moments: 

{Fa{t)) = 0, {a,b = x,y,z) 

and 

{Faih)Fb{t2)) = 2TjeSabSD{h - ^2), (2.12) 

where D = 2riQ represents the diffusion constant entering 
into a Fokker-Planck description. Choosing 8 to equal 
the initial energy then ensures that the average energy of 
the orbits remains unchanged. 

This equation is unsatisfactory here. Energy is con- 
served absolutely for frozen-A^ orbits, so that one must 
also impose energy conservation on any scheme which 
aims to mimic its effects. (For very small 77, energy re- 
mains almost conserved for very long times. However, 



comparatively small N should correspond to relatively 
large 77, which implies large changes in energy and, as 
such, significant changes in the phase space regions ac- 
cessible to the noisy orbit.) For this reason, the noisy 
integrations described here were performed using a mod- 
ified energy-conserving noise. 

What this entailed was (1) eliminating the dynamical 
friction altogether, (2) again imparting random kicks as 
in eq. (2.12), but (3) renormahsing the modified velocity 
at each time step by an overall factor, i.e., v(i -I- St) 
av{t + 5t), with a so chosen that E{t + 5t) = E{t). Mod- 
ulo this additional complication, the noise was integrated 
using a standard algorithm |Q based on a fourth order 
Runge-Kutta integration scheme with fixed time step 6t. 
Most of the integrations were performed for 5t = 2x 10~^, 
it having been confirmed that the statistical effects of 
the noise were insensitive to the precise value of 6t for 
St < 10~3. 



III. 'REGULAR' AND 'CHAOTIC ORBITS 

A. Ordinary Lyapunov exponents a useless 
diagnostic for macroscopic chaos 

In terms of a computation of Lyapunov exponents for 
individual orbits, it is difficult, if not impossible, to dis- 
tinguish between frozen- A'^ trajectories corresponding to 
regular orbits in a smooth potential and frozen- A'^ trajec- 
tories corresponding to smooth chaotic orbits. 'Regular' 
and 'chaotic' frozen- A^ orbits typically have positive Lya- 
punov exponents which are comparable in magnitude, 
and that magnitude is typically much larger than the 
magnitude of the largest Lyapunov exponent for chaotic 
orbits evolved in the smooth potential. 

This is illustrated by FIG. 1, the two panels of which 
exhibit estimates of the largest Lyapunov exponents for 
frozen- A^ realisations of the homogeneous ellipsoid - both 
with and without a central mass - and of the Dehnen 
mass distribution. For both distributions, estimates 
of X were computed for values of A'^ ranging between 
A^ = lO^-^ and A^ = 10^-^ Each curve in this FIGURE 
represents a mean value extracted by averaging over four 
different initial conditions. 

The three curves in the left panel correspond to regu- 
lar initial conditions evolved in an ellipsoid with M bh = 
(solid line) and 'sticky chaotic' (dashed lines) 

and 'wildly chaotic' (dot-dashed lines) initial conditions 
evolved in an ellipsoid with Mbh = lO""'^'^. Three points 
in this panel are apparent: (1) The estimates of x com- 
puted for these three different sets of initial conditions 
are very nearly equal in magnitude. (2) Consistent with 
the results described in Paper I, one observes no sys- 
tematic dependence on A^. (3) The typical value of x is 
comparatively large, much larger than the typical values 
Xs associated with motion in the smooth potential. For 
the values of energy E and mass AI bh used to generate 
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these x's, 'wildly chaotic' orbit segments in the smooth 
potential typically have xs 0.055 and 'sticky' segments 
have xs ~ 0.022. 

The facts that 'regular' orbits with Mbh = and 
'chaotic' orbits with Mbh > have x's that are ap- 
proximately equal, and that this value is much larger 
than xs, suggests strongly that these exponents reflect 
almost completely the effects of 'microscopic' chaos as- 
sociated with close encounters, and that the form of the 
bulk potential is largely immaterial. In both cases the 
frozen- A'^ orbits are moving through an ellipsoid with the 
same constant density; and, since the black hole mass 
Mbh — 10^^-^ ^ 1, the presence or absence of the black 
hole does not significantly impact the natural orbital time 
scale. 

The second panel in FIG. 1, generated for compara- 
tively low energies orbits evolved in the Dehnen poten- 
tial, shows an even closer agreement between initial 
conditions corresponding to regular and chaotic orbits. 
However this agreement does not directly corroborate the 
results derived from the first panel. The smooth Dehnen 
potential admits a coexistence of large measures of both 
regular and chaotic orbits and, as will be discussed in Sec- 
tion V, for particle number N < 10^, discreteness effects 
can trigger transitions between 'regular' and 'chaotic' be- 
haviour on time scales much shorter than the integra- 
tion times used to obtain estimates of the frozen-A'^ x- 
Over the time scales of interest it is impossible to make 
clean distinctions between 'regular' and 'chaotic' orbits 
in frozen-A^ realisations of the Dehnen distribution. At 
any rate, the frozen- A^ Lyapunov exponents are all much 
larger than the typical value xs ~ 0.1 for chaotic orbits 
evolved in the smooth Dehnen potential. 

B. The divergence of A^-body orbits from smooth 
characteristics 

One obvious way in which to effect macroscopic com- 
parisons of frozen-A^ orbits and smooth characteristics 
generated from the same initial condition is by computing 
such diagnostics as the configuration and velocity space 
separations, 

Dr{t)=\rs{t)^rN{t)\ and 
Dv{t)=\wsit)^VNit)l (3.1) 

where {rs,vs) and (rAr,vjv) denote, respectively, phase 
space coordinates for orbits in the smooth and frozen- A'^ 
density distributions. 

It was found in Paper I that, for individual regular 
orbits |p6|| , Dr and Dv typically grow linearly in time. 
A corresponding analysis for individual chaotic orbits 
is somewhat less conclusive since, apparently, different 
chaotic orbits can exhibit a larger degree of variety in 
their behaviour. However, by averaging over an ensem- 
ble of m different initial conditions, and computing the 
mean separation 



Vr{t) = -f2Dn, (3.2) 

i=l 

and an analogous 'Dv{t), one can again extract an unam- 
biguous trend. As for the case of regular orbits, so also for 
chaotic orbits, viewed macroscopically the quantities Vr 
and Vv diverge linearly in time. This is, e.g., illustrated 
in the first seven panels of FIG. 2, which was generated 
by effecting pointwise comparisons of 800 smooth and 
frozen- A^ orbits in the potential (2.6) for different values 
of N ranging between A^ = lO^-^ and A^ = 10^ ^ The 
initial conditions, chosen identically for each value of A^, 
were selected so as to correspond to 'wildly chaotic' or- 
bits, for which a typical short time Lyapunov exponent 
XS « 0.055. 

'Regular' and 'chaotic' frozen-A^ orbits are similar in 
the sense that they both diverge from a smooth charac- 
teristic with the same initial condition as a power law 
in time: Vr = A^t. However, the characteristic time 
scales involved are very different. For regular orbits in 
the Plummer potential, it was found (c/. eq. (4.8) in 
Paper I) that 

to ~ Ac^regN^'^tD, regular (3.3) 

where t^ a characteristic crossing time and Acreg is a 
constant of order unity. For chaotic orbits in the Dehnen 
potential and the ellipsoid plus black hole potential (2.6), 
one discovers instead that 

tG~AG,cha{\T).N)tD, chaotic (3.4) 

where, again, A^cha is of order unity. 

The goodness of fit to such a logarithmic dependence 
is exhibited in the final panel of FIG. 2, which exhibits 
growth times to derived from least squares fits to the 
data exhibited in the preceding panels. This panel should 
be contrasted with FIG. 8 in Paper I, which exhibits an 
analogous curve derived for regular orbits in the Plummer 
potential. 

C. Distinctions based on orbital 'complexity' 

The fact that, as A^ increases, frozen- A^ orbits re- 
main close to smooth characteristics with the same initial 
condition for progressively longer times implies that, in 
terms of visual appearance, they also tend to more closely 
resemble those smooth characteristics. This visual im- 
pression is easily corroborated by an examination of the 
complexity of the Fourier spectra, which characterise the 
extent to which the orbits are, or are not, nearly peri- 
odic. For comparatively small A^, frozen-A^ orbits corre- 
sponding to both regular and chaotic characteristic both 
look strongly aperiodic and 'wildly chaotic' in appear- 
ance. Not surprisingly, therefore, their complexities are 
very large, large compared even with the complexities 
associated with ordinary 'wildly chaotic' orbits evolved 
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in the corresponding smooth potential. However, as N 
increases, the complexities of regular and chaotic frozen- 
N orbits both decrease and, for sufficiently large N, the 
complexities n{N) converge towards the values ns asso- 
ciated with orbits in the smooth potential. For frozen- 
orbits corresponding to regular characteristics, this ng is 
typically quite small; for orbits corresponding to chaotic 
characteristics, ns is typically much larger. 

Analogous results obtain for orbits evolved in the pres- 
ence of ordinary friction and noise and, as such, it 
would seem natural to ask whether the observed varia- 
tions in complexity resulting from changes in N can also 
be mimicked by the energy-conserving noise considered 
in this paper. Overall, the answer to this would seem to 
be: yes. 

FIG. 3 exhibits the mean complexity n for represen- 
tative samples of orbits evolved in the potentials (2.2), 
(2.5), and (2.6). In each case, the complexities were com- 
puted from an orbital time series that sampled an orbit 
of duration T = 512 at intervals At — 0.05. The first two 
panels were generated for the Plummer and homogeneous 
ellipsoid potentials, both corresponding in the continuum 
limit to completely integrable motion. The last two pan- 
els exhibit complexities computed for initial conditions 
corresponding to wildly chaotic and 'sticky' chaotic orbits 
in the ellipsoid plus black hole potential. In each panel, 
the diamonds reflect complexities computed for frozen- 
N orbits with different values of N and the horizontal 
dashed line corresponds to the mean complexity for un- 
perturbed orbits with the same initial conditions evolved 
in the corresponding smooth potential. The triangles re- 
flect complexities computed for motion in the smooth po- 
tential perturbed by noise with = 1.0 and a coefficient 
rj related to N via a relation log;iQ t] = A — log^Q N. 

The constants A were not determined by using a least 
squares algorithm to make the two curves coincide, at 
least approximately. Rather, as described more care- 
fully in Section IV, the connection between r/ and N was 
effected by demanding that, for the case of regular or- 
bits, noise of amplitude 77 and discreteness effects asso- 
ciated with a system of N particles yield linear phase 
mixing at the same rates. For the Plummer potential, 
this requires Apium ~ 0.0. For the homogeneous ellip- 
soid, Aeiiip ~ 0.5. Because the black hole mass in the 
potential (2.6) is much smaller than the total mass of 
the system - Mbh = 10~^'^M - it should have only a 
minimal effect on the characteristic orbital time scales; 
and, for this reason, the last two panels involved fits as- 
suming A = AelUp ■ 

Overall, the correspondence between frozen-A^ and 
noisy complexities is quite good, except for very small N 
and very large 77 where systematic discrepancies can be 
seen to occur. What this means is that if the amplitude of 
the noise is tuned so as to reproduce the expected macro- 
scopic behaviour associated with regular phase mixing of 
orbit ensembles, individual noisy orbits - both regular 
and chaotic - will also manifest a degree of complexity 
comparable to that manifested by individual frozen-A^ 



orbits. 



IV. PHASE MIXING OF ORBIT ENSEMBLES 
A. Divergence of initially localised ensembles 

For particle number N < 10^ or so, initially localised 
'regular' and 'chaotic' ensembles both phase mix ex- 
tremely rapidly in a fashion that makes it virtually im- 
possible to distinguish between them. Inded, it would 
appear that, for such small N, mixing is dominated by 
discreteness effects and the form of the bulk potential 
is comparatively unimportant. However, for N > 10^'^ 
or so, it becomes possible to make clean distinctions be- 
tween the phase mixing of ensembles which correspond 
in the continuum limit to regular vis-a-vis chaotic orbits. 
It is, for example, evident that 'regular' ensembles dis- 
perse as a power law in time, whereas chaotic ensembles 
disperse in a fashion that is roughly exponential. 

Consider, e.g., the configuration dispersion ctj- associ- 
ated with an ensemble, which satisfies 

= (r^) - (r)^ (4.1) 

with 

{rn = -j:(^"+y"+^"r^"- (4-2) 

i—l 

For the case of a homogeneous ellipsoid with Mbh = 0, 
which corresponds in the continuum limit to integrable 
orbits all oscillating with the same natural frequencies, 
there is no phase mixing in the smooth potential, so 
that ar exhibits no systematic growth. If, however, the 
smooth potential is replaced by a frozen- A?^ potential, one 
discovers instead that ar grows as t^^"^. More precisely, 
the growth of the dispersion is well fit by a simple relation 
of the form 

ar = (^Ag)'/', (4.3) 

where 

taiN) = ANtD (4.4) 

and A is a constant of order unity. 

This ^^/^ behaviour is hardly surprising, corresponding 
as it does to the analytically predicted behaviour of an 
ensemble of 'noisy' orbits all evolved in a harmonic po- 
tential with the same natural frequencies (c/. 1^). The 
fact that to exhibits a roughly linear dependence on N 
is also expected, given the expectation |J] that discrete- 
ness effects should be manifested on a relaxation time 
tn cx N/ InN. The validity of the scaling ta ex: N is illus- 
trated in FIG. 7, which will be discussed more carefully 
below. 

That the growth in ar is slower than linear and that the 
growth time increases with increasing A^ are both evident 
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from FIG. 4, which exhibits cr^ for the same ensemble 
of initial conditions, evolved in frozen- A'^ potentials with 
variable particle number ranging between N = 10^'^ and 
N = 10^ 

For more generic models like the Plummer potential, 
where different integrable orbits oscillate with different 
frequencies, the behaviour is somewhat more complex: 
For smaller N the dispersion still grows as t^^^ on a time 
scale consistent with cq. (4.4), whereas, for larger N, 
the growth law is essentially linear. This behaviour can 
be interpreted as arising from a competition between two 
effects. Regular phase mixing, present even in the contin- 
uum limit, results in a dispersion that grows linearly in 
time; whereas an additional 'noise-induced phase mixing' 
typically induces a t^^^ divergence between perturbed 
and unperturbed orbits. For large A'', discreteness effects 
have only a small effect on the phase mixing, so that 
the frozen- iV ar exhibits the same linear dependence on 
time as the smooth potential ar- For smaller N discrete- 
ness effects become more important and induce a larger 
amplitude t^/^ growth. 

This behaviour is exhibited in FIG. 5, the six pan- 
els of which exhibit (Tr for evolution in the smooth 
Plummer potential and for frozen-A^ evolution with 
102-5 <N< 10^ It is evident that, for N > 10^ the 
overall growth is linear, at least until the dispersion 'sat- 
urates' at a value ar ^ 0.4, and that this linear evolution 
is nearly indistinguishable from the behaviour associated 
with the smooth evolution. The only obvious difference is 
that discreteness effects tend to 'fuzz out' the systematic 
oscillations associated with strictly periodic orbits in the 
unperturbed potential which yield the large 'spread' in 
values of ar superimposed on the overall linear growth. 
For smaller values of TV, ar{t) is better represented by 
a t^/^ growth law; and, quite apparently, the 'spread' in 
ar is substantially reduced. This reduction is an obvious 
manifestation of the fact that, for smaller N, the orbits 
in the ensemble exhibit significant deviations from peri- 
odicity, a fact manifested by the increased complexities 
discussed in the preceding Section. 

Implicit in the preceding is the assumption that dis- 
creteness effects really can be mimicked by noise. This 
was tested at two levels, namely (i) a qualitative, visual 
comparison of plots of ar{t) generated for both frozen- A?' 
and noisy ensembles, and (ii) a quantitative comparison 
of slopes associated with a = t/tc growth law. The 
degree to which noise can mimic discrete effects is evident 
visually from a comparison of individual panels in FIG. 5, 
which were generated from frozen- A'^ ensembles evolved 
in the Plummer potential, with the corresponding panels 
in FIG. 6, which were generated from the same ensembles 
of initial conditions, now evolved in the smooth potential 
in the presence of energy-conserving white noise. For par- 
ticle number as small as A^ = lO^ ^, the correspondence 
is comparatively poor, the frozen-A'' ar growing consid- 
erably more rapidly than the noisy ar- However, already 
for A^ = 10'^ the correspondence is quite reasonble and, 
for N > 10^, the frozen-A'' and noisy plots are virtually 



indistinguishable. 

The preceding comparison of frozen- A^ ensembles with 
fixed A^ and noisy ensembles evolved with 77 cx 1 /N was 
motivated by the fact that the noisy dispersions were also 
well fit by a t^^"^ growth law where, however, 

to = Bto/v, (4.5) 

with B a constant of order unity. A comparison of eqs. 
(4.4) and (4.5) implies a natural identification 

logio v = A- logio N, (4.6) 

with A yet another constant. For the ellipsoid potential, 
the best tit A = O.OibO.l. For the Plummer potential 
A = 0.6±0.1. The identification between FIGS. 5 and 
6 was effected assuming A = 0.5. The extent to which 
the growth rates for frozen-A^ and noisy orbits can be 
related by the simple relation (4.6) can be gauged from 
FIGS. 7, which superimpose plots of tciN) and toiri) 
for the Plummer and ellispoid potentials, with N and r] 
related by eq. (4.6). 

The behaviour exhibited by ensembles corresponding 
to chaotic orbits is very different. In this case, the disper- 
sion associated with an ensemble evolved in the smooth 
potential typically diverges exponentially in time, and 
discreteness effects only serve to accelerate this growth. 
This is, e.g., evident from FIG. 8, which plots \nar for 
an ensemble of initial conditions evolved in the potential 
(2.6) that corresponds in the continuum limit to 'wildly 
chaotic' orbits. It is apparent that the growth of ar is 
not strictly exponential, but it certainly is faster than 
the exponential growth associated with the smooth po- 
tential, which is exhibited in the final panel. One other 
point is also evident: For the case of the Plummer poten- 
tial, a frozen- A^ integration with N as small as A^ ~ 10^'^ 
yields phase mixing almost identical to that observed in 
the smooth potential. For the case of the chaotic ensem- 
ble evolved in the potential (2,6), even A'' = 10^-^ yielded 
phase mixing that was much more rapid than that asso- 
ciated with the corresponding smooth potential. 

FIG. 10 exhibits an analogous plot, again generated 
for the potential (2.6), but now considering an ensemble 
of initial conditions which, in the continuum limit, cor- 
respond to comparatively 'sticky' chaotic orbits which 
initially disperse much more slowly. (The typical value 
of the largest short time Lyapunov exponent for smooth 
orbits in this 'sticky' ensemble was x(t = 256) ~ 0.22. 
The typical value for the 'wildly chaotic' ensemble was 
X ~ 0.55.) For A^ < 10^ or so, FIGS. 8 and 10 are com- 
paratively similar, the 'stickiness' manifested in the con- 
tinuum limit being largely lost. However, for larger N 
more conspicuous differences become apparent. It is, 
for example, clear that, for particle number as large as 
A'^ = 10^'^, the dispersion for the 'sticky' ensemble only 
becomes 'macroscopic' on a time scale appreciably longer 
than the time scale for the 'wildly chaotic' ensemble. 

To a considerable degree, the behaviour exhibited by 
chaotic frozen-A/' ensembles is again well-mimicked by 
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energy-conserving noise. This is, e.g.. evident from 
FIGS. 9 and 11 which, respectively, exhibit the same 
ensembles of initial conditions as FIGS. 8 and 10, now 
evolved in the smooth potential in the presence of energy- 
conserving white noise. The first seven panels in each of 
these FIGURES involved the same identification between 
N and r) as did FIG. 6, this motivated by the recog- 
nition that the black hole mass Mbh = 10~^'^ is too 
small to significantly alter the natural dynamical time 
to- The final panel in each FIGURE exhibits Ina^ for 
T] = 10~^'^, the largest value of rj that does not result in 
phase mixing that is significantly more rapid than that 
associated with motion in the unperturbed smooth po- 
tential. Presuming that the correspondence between TV 
and 7] established here can be extrapolated to larger A'' 
and smaller r], it follows that, for the ensembles consid- 
ered here, one would require particle number as large as 
N ~ 10^ before discreteness effects become unimportant 
over the time scales of interest! 

One obvious issue is the extent to which the behaviour 
of orbit ensembles evolved in different frozen-A^ realisa- 
tions of the same potential with the same number N is, 
or is not, the same. For the case of the intcgrable Plum- 
mer and ellipsoid potentials, it would appear that, even 
for numbers as small as A/' ^ 10'^, the statistical proper- 
ties of orbit ensembles are essentially the same for differ- 
ent frozen- A'' realisations. This, however, is not true for 
chaotic ensembles. In this case, one appears to require 
A'^ > 10^ or so before noticeable differences between dif- 
ferent realisations are suppressed. This is, e.g., evident 
from FIG. 12, which exhibits cr^ for the same initial con- 
ditions used to generate FIGS. 10 and 11, now focusing 
on a time interval only one fifth as long and including 
with dotted lines the results of the different realisations 
which were averaged to yield the solid line. 

B. Approach towards a 'well-mixed' state 

The fact that initially localised ensembles of chaotic 
orbits evolved in a smooth potential should diverge ex- 
ponentially is more or less obvious; and the simulations 
described in this paper show that graininess associated 
with finite A'' has the same qualitative effect as noise 
with characteristic amplitude rj oc 1/N. Frozen- A'^ evo- 
lution results in a divergence which is even faster than 
that associated with the unperturbed ensemble, albeit no 
longer strictly exponential. 

Less obvious, but also true, is the fact that chaotic en- 
sembles integrated in a smooth potential tend to evolve 
towards a near-equilibrium state, and that this evolution 
typically proceeds exponentially in time. For systems 
manifesting the reflection symmetries associated with the 
potentials explored in this paper, this implies in partic- 
ular that moments like (x) and (vx) tend to zero ex- 
ponentially. The obvious question then is whether dis- 
creteness effects, which can significantly accelerate the 



initial rate at which ensembles disperse, also act to ac- 
celerate the rate at which ensembles evolve towards a 
near-equilibrium. 

Similarly, one can consider the effects of graininess on 
regular orbits. Unperturbed regular orbits evolved in a 
smooth potential do not exhibit a rapid exponential ap- 
proach towards a (near-)equilibrium corresponding to a 
finite phase space volume. Rather, what one observes 
is a more modest power law evolution towards a near- 
uniform population of the invariant tori to which they 
are restricted. Allowing for discreteness effects allows the 
orbits to escape from these invariant tori and, as such, 
one might again ask: is there an approach towards some 
more general 'near-equilibrium' ? And, if so, how does 
this approach proceed in time? 

For the case of potentials which, in the continuum 
limit, correspond to integrable systems, one finds that 
'weak' discreteness effects associated with comparative 
large N have only a minimal effect on such moments 
as (x). If, however, one allows for smaller A^ systems, 
where discreteness effects become more important, these 
moments will in fact exhibit an unambiguous exponen- 
tial convergence towards zero. This is, e.g., illustrated 
in FIGS. 13 and 14, which exhibit the quantity ln|(z)| 
for frozen- A^ ensembles evolved, respectively, in the inte- 
grable ellipsoid and Plummer potentials. The fact that 
this convergence seems to terminate at In | (2) | ~ — 3 is 
a finite sampling effect, reflecting the fact that the en- 
sembles were comprised of only 1600 orbits. Even if one 
were to select 'at random' 1600 points from a continuous 
distribution with {z) = 0, one would have a sample for 
which {z) ^ 0. 

For the case of chaotic potentials one observes an ex- 
ponential decrease in such moments, even in the con- 
tinuum limit, but discreteness effects again serve to in- 
crease the rate associated with this evolution. This is, 
e.g., illustrated in FIGS. 15 and 17, which exhibit In 
for the 'wildly chaotic' and 'sticky' ensembles used to 
generate FIGS. (3.5) and (3.7). The degree to which 
this enhanced exponential evolution can be mimicked by 
energy-conserving noise is illustrated in FIGS. 16 and 
18. As for FIGS. (3.6) and (3.8), the final panel in each 
of these FIGURES corresponds to noise with amplitude 
7] = 10~^'^ which, once again, is the weakest noise not 
to occasion a significant increase in the overall eflacacy of 
chaotic mixing. 



V. TRANSITIONS BETWEEN 'REGULAR' AND 
'CHAOTIC BEHAVIOUR 

It is well known that friction and noise can convert a 
regular orbit into a chaotic orbit and vice versa. Suppose, 
e.g., that an initially chaotic orbit is evolved in a smooth 
potential in the presence of noise for some finite period 
and that the noise is then turned off. If the evolution is 
continued in the absence of noise, one may then find that 



9 



the orbit has become regular with no positive Lyapunov 
exponents. Why this can happen is easy to understand: 
Noise serves to continually 'bump' the orbit from one 
smooth characteristic to another and it is quite possible 
that such 'bumps' will eventually deflect the orbit from 
a chaotic to a regular characteristic. To the extent that 
discreteness effects can be mimicked by noise, one would 
anticipate the possibility of similar transitions in an iV- 
body evolution p9| . 

This was not an issue for orbits evolved in the po- 
tentials (2.2) and (2.5), both of which are integrable. 
This possibility is also unimportant for the nonintegrable 
black hole plus ellipsoid potential discussed above since, 
for the energies and black hole masses that were consid- 
ered, the energetically accessible phase space is almost 
completely chaotic. However, the possibility of such tran- 
sitions is a major issue for the Dehnen potential where, 
for most energies, the constant energy hypersurface ad- 
mits large measures of both regular and chaotic orbits. 

One hint that such transitions might be present derives 
from a computation of I?r, the mean separation between 
frozen- iV orbits and smooth characteristics with the same 
initial condition, generated from ensembles of initial con- 
ditions corresponding to only regular or only chaotic or- 
bits. One example is provided in FIG. 19, which exhibits 
Dr{t) for ensembles evolved in frozen-A^ realisations of 
the Dehnen distribution with iV = 10"*, iV = lO'*-^, and 
N = 10^. The three left panels were generated from regu- 
lar initial conditions; the three right panels from chaotic 
initial conditions. Both sets of initial conditions were 
selected to have very low energies, so that they were re- 
stricted to the central portions of the system [g8[ . 

It is clear that the chaotic initial conditions yield a 
more 'irregular' time dependence although, as expected, 
in both cases Dr evidences a nearly linear growth in time. 
However, it is also obvious that the growth time to for 
the chaotic initial conditions is not much longer than the 
growth time for the regular initial conditions. For all 
three values of N, the best fit value of to is comparable 
for the regular and chaotic initial conditions although, 
for N = 10^, it is clear that Dr grows somewhat more 
slowly for the regular orbits. The obvious interpretation 
is that, even at very early times, a significant fraction 
of the initial conditions have switched between 'regular' 
and 'chaotic' behaviour. 

Another, even more striking, indication that such tran- 
sitions occur is provided by computing distributions of 
orbital complexities. This is illustrated in FIG. 20, which 
exhibits the complexity n for the low energy ensembles 
used to generate FIG. 19, as well as for a comparable 
set of integrations performed for ensembles with substan- 
tially higher energies The three left panels corre- 
spond to the lower energy initial conditions; the three 
right panels to the higher energy initial conditions. In 
each case, the chaotic orbits are represented by a solid 
curve and the regular orbits by a dashed curve. 

It is clear that, for the lower energy evolutions, the 
distributions of complexities for the regular and chaotic 



initial conditions are almost identical, although there is 
some hint that, for N = 10^, the regular orbits are some- 
what 'less complex.' However, for the higher energy evo- 
lutions regular and chaotic initial conditions yield very 
different distributions of complexity. For all three val- 
ues of N, the distribution for the chaotic ensembles is 
unimodal and peaked at a value n ^ 500 — 600. For the 
two larger values of N, the distributions for the regu- 
lar orbits are also singly peaked, but the peaks occur at 
significantly lower values n ^ 200. Significantly, though, 
for the lowest value of N the distribution of complexities 
for the regular initial conditions is distinctly bimodal, 
with one peak at a value similar to that obtained for the 
chaotic initial conditions and another at a significantly 
smaller value. This sort of bimodal distribution is ex- 
actly what one would expect for an ensemble of orbits 
which is partially regular and partially chaotic [^7|] . 

Direct proof that transitions between regularity and 
chaos occur, and an estimate as to the time scale on 
which such transitions occur, is straightforward to ob- 
tain. Given orbital data from frozen-A'' ensembles in- 
tegrated for a total time T, one can use snapshots at 
earlier times t < T to generate new ensembles of initial 
conditions and integrate those ensembles in the smooth 
potential to compute Lyapunov exponents. Comparing 
the results of such integrations with integrations of the 
original initial conditions in the smooth potential then 
permits a concrete probe of the extent to which discrete- 
ness effects have made a regular ensemble 'less regular' 



and/or a chaotic ensemble 'less chaotic' |30 



The results of such an analysis are summarised in 
FIG. 22, which plots the mean short time Lyapunov ex- 
ponent (xs) for smooth orbit ensembles as a function 
of T. The three left panels were generated for the low 
energy ensembles; the three right panels for the high en- 
ergy panels. In each case the initially chaotic ensembles 
are represented by diamonds and a solid line; the ini- 
tially regular ensembles are represented by triangles and 
a dashed line. The three rows again correspond respec- 
tively to iV = lO'^, N = 10''■^ and N = 10^. 

In each case it is clear that, overall the mean xs for the 
initially chaotic ensembles decreases systematically with 
T and that the mean xs for the initially regular ensem- 
bles increases. This is of course exactly what would be 
expected if, as r increases, progressively larger numbers 
of transitions between regularity and chaos have occured 
during the frozen- evolution. Indeed, to the extent that 
the two ensembles of initial conditions sample the same 
phase space regions, one would expect that the values of 
Xs for the regular and chaotic ensembles should converge 
to a common value. It is evident that, for the lower en- 
ergy ensemble, convergence or near-convergence has in 
fact been achieved for the largest values of r, and it is 
also apparent that convergence happens more rapidly for 
smaller N. This again is exactly as expected. Larger N 
corresponds to weaker 'noise,' but weaker noise should 
be less effective in inducing transitions. 

A priori it might seem surprising that the mean xs 
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for the regular orbits approaches its limiting value much 
more quickly than does the mean for the chaotic orbits. 
This, however, is likely a manifestation of the phase space 
structure associated with the Dehnen potential. For both 
low and high energies, the initially regular ensemble cor- 
responds in the smooth potential to box orbits which 
occasionally pass quite close to the center of the system; 
but, in the center, the regular and chaotic phase space 
regions are entangled in a very complex fashion, so that 
it is comparatively simple for initially regular orbits to 
be deflected to chaotic trajectories. 

VI. DISCUSSION 

The experiments described in this paper lead to several 
unambiguous conclusions regarding the behaviour of or- 
bits and orbit ensembles evolved in frozen- AT realisations 
of smooth density distributions. 

It is, for example, clear that ordinary Lyapunov ex- 
ponents computed for individual frozen-iV trajectories 
do not provide a useful characterisation of the degree 
of macroscopic chaos manifested by these trajectories. 
Different initial conditions with the same energy evolved 
in frozen- A'' realisations of a specified potential typically 
have Lyapunov exponents that are comparable in mag- 
nitude and exhibit little if any dependence on N, even 
though, in terms of their bulk properties, one may look 
very nearly regular and another wildly chaotic. 

This can be interpreted by asserting that, at least for 
large N, one can make comparatively clear distinctions 
between two types of chaos which may be associated with 
the AT-body problem. On the one hand, there is micro- 
scopic chaos associated with close encounters between 
nearby masses. This chaos, which is presumably respon- 
sible for the large Lyapunov exponents associated with 
frozen- A'^ orbits, is (almost) always present but, at least 
for comparatively large A^, tends to 'wash out' macro- 
scopically. On the other hand, there is the possibility 
of macroscopic chaos which, if present in the continuum 
limit, will also have manifestations in frozen- A^ simula- 
tions. 

The experiments also demonstrate a clear sense in 
which, as N increases, frozen-A?^ trajectories become 
'more similar' to smooth potential characteristics gen- 
erate from the same initial condition. In particular, this 
implies that it becomes progressively easier to distinguish 
between 'regular' and 'chaotic' macroscopic behaviour. 
This similarity can be quantified in at least three differ- 
ent ways: 

1. The rate at which frozcn-N orbits and smooth char- 
acteristics diverge. For both regular and chaotic initial 
conditions, frozen- A^ trajectories and smooth character- 
istics tend to diverge linearly in time. However, the 
A/'-dependence of the time scale to associated with this 
divergence differs dramatically for regular and chaotic 
orbits. For initial conditions corresponding to regular 



orbits, to ex N^^^. For initial conditions corresponding 
to chaotic orbits, tc oc InA^. It follows that, for large 
N, frozen- A'' trajectories and smooth characteristics cor- 
responding to regular orbits remain 'close' much longer 
than do trajectories and characteristics corresponding to 
chaotic orbits. 

2. The 'complexity' of Fourier spectra constructed from 
orbital time series. For small A'', both regular and chaotic 
initial conditions will, when integrated into the future, 
yield Fourier spectra that are much more complex than 
the spectra associated with an evolution in the smooth 
potential. In particular, both regular and chaotic orbits 
will yield spectra of comparable complexity. However, as 
A^ increases the complexities decrease and, for sufficiently 
large A^, one sees a convergence towards the complexities 
appropriate for smooth characteristics with the same ini- 
tial conditions, be these either regular or chaotic. 

3. The bulk properties of phase mixing, as probed, e.g. by 
lower order moments. Phase mixing of initially localised 
ensembles in frozen- A^ systems is invariably more efficient 
than phase mixing in the corresponding smooth poten- 
tial but, as N increases, the observed evolution comes to 
more closely resemble phase mixing of the same initial 
conditions evolved in the smooth potential. 

At least for comparatively large N, many discreteness 
effects can be well mimicked by energy-conserving white 
noise with amplitude rj oc 1/N. This is in close agree- 
ment with naive expectations based on the modeling of 
discreteness effects as a sequence of incoherent binary en- 
counters, which would suggest q oc (In N)/N; and indeed, 
the simulations arc also consistent with the latter depen- 
dence. Modeling discreteness effects in terms of noise 
works well both in predicting the expected complexities 
of individual frozen- A^ trajectories and moments associ- 
ated with phase mixing. This indicates that, at least 
for large A'^ noise can be used to model both bulk statis- 
tical properties of orbit ensembles and qualitative prop- 
erties of individual orbits. This suggests strongly that 
investigations of the effects of white noise on orbits in 
nonintcgrable potentials - which are much less expensive 
computationally than frozen- AT integrations - can be in- 
terpreted as providing information about the effects of 
graininess. 

Nevertheless, it is clear that modeling discreteness ef- 
fects by white noise is not completely satisfactory. The 
agreement between frozen- A^ and noisy integrations is 
poor for comparatively small A^; and noisy integrations 
of the form described in this paper cannot be used to 
obtain estimates of the Lyapunov exponents associated 
with frozen- A^ orbits. The development of a 'more real- 
istic' noise and its use to model systems with N ^ 10^ 
is currently underway. 

Discreteness effects can induce transitions between 
different orbit 'types,' including both transitions be- 
tween regular and chaotic behaviour, which are impos- 
sible in the smooth potential, and transitions between 
(say) 'sticky' and 'wildly chaotic' behaviour which, in 
the smooth potential, is not forbidden but typically ce- 
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curs only on a much longer time scale. Not surprisingly 
these transitions appear to be more common for smaller 
N; and indeed, the simulations arc consistent with the 
interpretation that, for sufficiently large A'', transitions 
between regular and chaos become essentially impossi- 
ble and transitions between sticky and wildly chaotic be- 
haviour happen no more quickly than in the continuum 
limit. 

Finally, and perhaps most importantly, it would ap- 
pear that, at least in terms of macroscopic properties, it 
does make sense to speak of a smooth N ^ oo continuum 
limit. However, convergence towards this limit is much 
slower for density distributions which, in the continuum 
limit, correspond to nonintegrable potentials that admit 
chaotic orbits. 
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FIG. 1. (a) Estimates of the largest Lyapunov exponent 
for orbits evolved in frozen-A'^ realisations of a homoge- 
neous ellipsoid: integrable initial conditions evolved with 
Mbh ~ (solid line), 'sticky' initial conditions evolved with 
AIbh ~ 10"*'"" (dashed line), and 'wildly chaotic' initial con- 
ditions with Mbh = lO"'^'^ (dot-dashed line), (b) Lyapunov 
exponents evolved in frozen-A realisations of the Dehnen po- 
tential for low energy regular (solid line) and chaotic (dashed 
line) initial conditions. 




FIG. 2. The mean separation, Vr, between frozen-A or- 
bits and smooth characteristics with the same initial con- 
ditions, computed for ensembles of 800 chaotic initial con- 
ditions evolved in the potential (2.6) for variable A. (a) 
A = 10^'^ (b) A = 10^ (c) A = 10='■^ (d) A = 10'*. 
(e) A = 10*'^ (f) A = 10^ (g) A = 10^'^ (h) The growth 
rate to (A) extracted from the preceding panels, assuming a 
linear growth law. The solid line overlays a least squares fit 
tG = A^B logio A. 
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FIG. 3. The mean complexity n for orbits evolved in the 
potentials (2.2), (2.5), and (2.6), considering both frozen-iV 
orbits with variable N (diamonds) and smooth orbits per- 
turbed by noise with G = 1.0 and variable rj (triangles), rj 
was related to A'' by the relation r] = e^/N, with A deter- 
mined as described in the text. The dashed horizontal line 
exhibits the mean complexity of orbits with the same initial 
conditions evolved in the smooth potential in the absence of 
noise, (a) Regular orbits in the Plummer potential, (b) Reg- 
ular orbits in the homogeneous ellipsoid potential, (c) Wildly 
chaotic orbits in the ellipsoid plus black hole potential, (d) 
'Sticky' chaotic orbits in the ellipsoid plus black hole poten- 
tial. 
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FIG. 4. The configuration dispersion ar associated with an 
initially localised ensemble evolved in frozcn-A^ realisations 
of the integrable ellipsoid potential (2.5) for variable A'", (a) 
N = 10^-^ (b) N = 10^. (c) N = 10^-^ (d) N = 10*. (e) 
N = 10" ^ (f) N = IQ^. 
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FIG. 5. The configuration dispersion Ur associated with an 
initially localised ensemble evolved in frozen- A'' realisations of 
the integrable Plummer potential (2.2) for variable N . (a) 
iV = 10^ (b) N = 10^ ^ (c) N = 10*. (d) N = 10* ^ (e) 
N = 10'' . (f) Evolution in the smooth potential. 
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FIG. 6. The configuration dispersion Cr associated with an 
initially localised ensemble evolved in the smooth Plummer 
potential (2.2), perturbed by noise with Q = 1.0 and vari- 
able 7?. (a) T] = 10"^ (b) 77 = 10-^-^ (c) T] = 10"*. (d) 
ri = 10-*-^ (e) r? = 10"^ (f) r] = 10"^ ^ 
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FIG. 7. (a) The growth time tG{N) for the dispersion cTr for 
orbit ensembles evolved in frozen- A?^ realisations of the Plum- 

mer potential (diamonds) and the corresponding growth time 
tG{ri) for orbit ensembles evolved in the smooth Plummer po- 
tential in the presence of energy-conserving white noise with 
O = 1.0 and logj^,, — — log^gN (triangles), (b) tG{N) and 
tairi) for orbit ensembles in the homogeneous ellipsoid poten- 
tial, now setting logj^Q r] = — logj^Q AT -|- 0.6. 
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FIG. 8. The configuration dispersion ffr associated with an 

initially localised ensemble of wildly chaotic orbits evolved in 
frozen- TV realisations of the nonintegrable potential (2.6) for 
variable N. (a) = 10^ ^ (b) N 
N = 10^. (e) N = W^-^. (f) N = 
Evolution in the smooth potential. 
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FIG. 9. The configuration dispersion Ur associated with an 
initially localised ensemble of wildly chaotic orbits evolved 
in the nonintegrable potential (2.6) but perturbed by en- 
ergy-conserving white noise with G = 1.0 and variable rj. (a) 
r; = 10 
r; = 10- 



(b) r; = 10-^ ^ (c) r) = 10"^ (d) r? = 10-=^ ^ (e) 

(f) n = 10-' ^ (g) V = lo-^ (h) V = w-'-'- 
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FIG. 10. The configuration dispersion (Jr associated with 
an initially localised ensemble of 'sticky' chaotic orbits evolved 
in frozen-iV realisations of the nonintegrable potential (2.6) 
for variable N. (a) N = 10^ ^ (b) AT = 10^ (c) N = 10^■^ 
(d) N = 10"*. (e) iV = lO"* ^ (f) N = 10^ (g) iV = 10^ ^ (h) 
Evolution in the smooth potential. 
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FIG. 11. The configuration dispersion Ur associated with 
an initially localised ensemble of 'sticky' chaotic orbits evolved 
in the nonintegrable potential (2.6) but perturbed by noise 
with e = 1.0 and variable 77. (a) r; = 10"^. (b) 77 = lO"^'^. 
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(e) 7? = 10-''. (f) 77 = 10" 
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FIG. 12. The configuration dispersion ar associated with 
the ensemble exhibited in the preceding FIGURE, now re- 
stricted to a shorter time interval, (a) = 10^'^. (b) 
N = 10^-^ (c) N = 10* ^ (d) N = lO^-^ 
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FIG. 13. The quantity ln|(2:)| for an initially localised en- 
semble evolved in frozen- A'^ realisations of the integrable el- 
lipsoid potential (2.5) for variable N. (a) N = 10^ ^ (b) 
N = W^. (c) = 10^■^ (d) AT = 10". (e) = 10'' ^ (f) 
AT = 10^ 



FIG. 14. The quantity lnl(z)| for an initially localised 
ensemble evolved in frozen-A^ realisations of the integrable 
Plummer potential (2.2) for variable N. (a) A^ = 1000. (b) 
N = W^-\ (c) N = 10". (d) N = 10"■^ (e) N = 10^ (f) 
Evolution in the smooth potential. 
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FIG. 15. The quantity In for an initially localised en- 
semble of wildly chaotic orbits evolved in frozen-A'^ realisa- 
tions of the nonintegrable potential (2.6) for variable N . (a) 
N = 10^■^ (b) N = 10^ (c) iV = 10^-'''. (d) = 10*. (c) 
N = 10''■^ (f) N = 10^. (g) N = 10^ '\ (g) Evolution in the 
smooth potential. 
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FIG. 16. The quantity ln|(2;)| for an initially localised en- 
semble of wildly chaotic orbits evolved in the nonintegrable 
potential (2.6) but perturbed by noise with G = 1.0 and vari- 
able r?. (a) 7j = 10~^ (b) 77 = 10"^ ^ (c) ?? = 10"'\ (d) 
V = W-''-'. (e) n = W-\ (f) V = 10-*^ (g) V = 10-^ (h) 
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FIG. 17. The quantity ln[(2;)| for an initially localised en- 
semble of 'sticky' chaotic orbits evolved in frozen- A'^ realisa- 
tions of the nonintegrable potential (2.6) for variable N. (a) 
N = 10^■^ (b) N = 10^ (c) iV = 10^-'''. (d) = 10*. (c) 
N = 10''■^ (f) N = 10^. (g) N = 10^ '\ (g) Evolution in the 
smooth potential. 
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FIG. 18. The quantity ln|(«)| for an initially localised en- 
semble of 'sticky' chaotic orbits evolved in the nonintegrable 
potential (2.6) but perturbed by noise with G = 1.0 and vari- 
able r?. (a) 7j = 10~^ (b) 77 = 10"^ ^ (c) ?? = 10"'\ (d) 
V = W-''-'. (e) n = W-\ (f) V = 10-*^ (g) V = 10-^ (h) 




FIG. 19. The quantity Vr for two different initially lo- 
calised ensembles of 'sticky' chaotic orbits with the same (low) 
energies, each evolved in frozen- iV realisations of the Dehnen 
potential, (a) A 'regular' ensemble with N — 10''. (b) A 
'chaotic' ensemble with — 10*. (c) The 'regular' ensemble 
with N = 10" ". (d) The 'chaotic' ensemble with N = lO"*'^. 
(a) The 'regular' ensemble with N = W^. (b) The 'chaotic' 
ensemble with AT = 10^. 
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FIG. 21. The mean Lyapunov exponent (xs), for ensem- 
bles evolved in the smooth Dehnen potential, selecting as ini- 
tial conditions the final phase space coordinates of orbits that 
had been evolved in frozen-TV potentials for some time r. (a) 
Initially regular (solid line) and chaotic (dashed line) low en- 
ergy orbits evolved with N = W^. (b) Initially regular (solid 
line) and chaotic (dashed line) higher energy orbits evolved 



with N = 10*. (c) The same as (a) but for N 
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FIG. 20. The distribution of complexities, N{n), for two 
different initially localised ensembles of 'sticky' chaotic orbits 
with the same (low) energies, each evolved in frozen-A^ reali- 
sations of the Dehnen potential, (a) A 'regular' ensemble with 
N = 10". (b) A 'chaotic' ensemble with A^ = lO". (c) The 
'regular' ensemble with A^ = lO"'^. (d) The 'chaotic' ensem- 
ble with AT = 10"-^. (e) The 'regular' ensemble with AT = 10^. 
(f) The 'chaotic' ensemble with AT = 10^. 



